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Multiple  Nonlinear  Diffusion  Equation 

In  the  previous  report  we  proposed  Nonlinear  Diffusion  Equation  ( NDE ) 
which  is  diffusion  approximation  to  Hasselmann  equation  (HE)  for  gravity 
waves.  NDE  is  the  second-order  partial  differential  equation  describing  the 
diffusion  of  the  density  of  action  for  gravity  waves  n(k,  <f>) 

Tt=TTkkh Tkk'V  +  Wn*  +  Tn’n  =  n r  =  r<*' «  W 

where  k  and  <f>  are  the  modulus  of  wavenumber  and  polar  angle  in  Fourier 
space,  T(k)  is  the  coefficient  responsible  for  external  forcing  and  viscous 
damping;  a  is  the  constant  to  be  found  from  the  comparison  of  numerical 
simulation  of  NDE  with  numerical  simulation  of  HE. 

In  the  absence  of  external  forcing  and  viscosity  NDE  has  the  same  mo¬ 
tion  integrals  as  HE  -  wave  action,  energy  and  momentum;  it  also  has  similar 
to  HE  asymptotic  Kolmogorov  spectra  corresponding  to  given  energy,  mo¬ 
mentum  and  wave  fluxes. 

After  change  of  variables  u>  =  k*  equation  (1)  takes  the  form 


dn 

dt 


1  d2  1  d2  ' 
2duj2  +  u)2d<j)2  _ 


u2in3  +  Tn 


The  procedure  of  comparison  of  NDE  and  HE  consists  in  comparison  of  the 
right-hand  side  (collision  term)  for  two  model  being  calculated  on  JONSWAP 
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spectrum.  Results  of  such  calculation  for  the  equation  (2)  have  shown  an  ex¬ 
cellent  agreement  with  the  results  of  HE  in  directional  distribution.  Agree¬ 
ment  between  angular  distributions  was  not,  however,  very  good  [2]. 

To  improve  an  agreement  between  HE  and  NDE  we  propose  Multiple 
Nonlinear  Diffusion  Equation  ( MNDE ) 


dn 

~di 


a 


Id2  Id2 

+ 


2  div2  uj 2  d(j)2 


u24(ain3  +  a2n2(n)  -1-  a3n(n2)  +  0'4n(n)2)  +  Tn  (3) 


where  brackets  mean  angle  averaging  in  ^-direction: 


1  r** 

{n)  =  2ii 

This  equation  is  identical  to  the  equation  (2)  if  04  ^  0  and  a*  =  0,  i  =  2, 3, 4. 
Extra  diffusion  terms  model  non-locality  of  interaction  in  HE  (that’s  where 
the  word  “multiple”  stem  from  in  MNDE). 

The  value  of  the  constants  oii  (i  =  1,2, 3, 4),  as  usual,  has  to  be  defined 
from  the  comparison  of  the  results  of  calculation  of  the  right-hand  sides  for 
MNDE  (3)  and  HE  being  calculated  on  the  JONSWAP  spectrum.  Such 
multi-parameter  optimization  procedure  is  nontrivial  and  is  the  subject  of 
future  research. 

In  the  current  report  we  represent  the  results  of  “informal”  choice  of  the 
coefficients  a*  (i  =  1,2, 3, 4)  which  demonstrate  that  even  simple  modifica¬ 
tion  of  the  model  (2)  gives  dramatically  better  agreement  in  the  angular 
distribution  between  HE  and  NDE. 

We  developed  numerical  solver  of  the  equation  (3)  allowing  to  calculate 
the  temporal  evolution  of  the  gravity  surface  wave  spectra.  To  get  the  “feel¬ 
ing”  of  the  effect  of  the  terms  in  the  right-hand  side  of  (3)  corresponding  to 
different  aij,  (z  =  1,2,3, 4)  we  have  made  four  independent  runs  of  MNDE 
solver  corresponding  to  the  sets: 

begin  eqnarray*  a,\  =  1,  a2  —  0,  a3  =  0, 04  =  0 

Q!i  =  0,  a2  =  1,  CK3  =  0, 0:4  =  0 
ai  =  0,  a2  —  0,  a3  =  1,  aq  =  0 
ai  —  0,  a2  =  0,  a3  =  0, 04  =  1 

It  is  interesting  that  three  first  runs  gave  close  results  (up  to  the  scaling 
constant) .  This  fact  prompted  us  to  make  the  conclusion  that  simple  choice 
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of  cti  ^  0, 0.2  =  0,  a3  =  0, 04  /  0  could  produce  better  agreement  in  angular 
distribution  with  the  results  of  simulation  of  H E  by  Resio  and  Tracy  [2]. 
The  choice  of  the  constants  ot\  =  0.467- 10-2  and  a4  =  0.14  gives  directional 
and  angular  distributions  represented  on  the  Fig.l  and  Fig.2  (solid  line  cor¬ 
responds  to  the  results  of  MNDE  and  crosses  to  the  results  of  HE  by  Resio 
and  Tracy) .  It  is  seen  that  an  agreement  is  excellent  for  directional  spectrum 
and  quite  good  for  angular  distribution  at  /  =  0.148. 

We  hope  to  make  the  angular  distribution  correspondence  even  better  as  a 
result  of  calculation  of  the  coefficients  ct;  ( i  =  1, 2, 3, 4)  from  multi-parameter 
optimization. 
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